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Separation of dissipation from diffusion 

Tong Zhou 

Department of Physics, University of Cahfornia, Santa Barbara, CA 93106 
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We study velocity correlations induced by diffusion and dissipation in a simple dissipative dynam- 
ical system. We observe that diffusion, as a result of time reversible microscopic processes, leads 
to correlations with different spatial parity from those caused by dissipation, consisting of time 
irreversible microscopic processes. The velocity correlations observed conflict with the "molecular 
chaos" assumption of the Boltzmann equation. To account for the apparent correlation structures 
in our simple model system, we propose that they indicate geometric distortion of the phase space. 
PACS numbers: 05.70.Ln 05.60. Cd 02.40.-k 



The second law of thermodynamics states that a sys- 
tem left to itself will evolve towards a state with maxi- 
mum entropy, i.e. maximum randomness. We are famil- 
iar with order due to interaction where, for example for 
systems in equilibrium at certain temperature, different 
strength of interaction among microscopic constituents 
leads to states with different degrees of order — solid, liq- 
uid, or gas. In this letter we study a different type of 
order, order purely induced by diffusion and dissipation. 

Nonequilibrium phenomena have long attracted much 
research interest. The Boltzmann equation is probably 
the most celebrated theory, yet its crucial assumption of 
"molecular chaos" needs scrutiny before it can be ap- 
plied to specific systems, see for example Linear 
response theory also has wide applications in con- 
densed matter physics. Some of recent studies include 
derivation of Ohm's Law for the periodic Lorentz gas 
the connection between the smoothness of distribution 
function and the number of degrees of freedom fi- 
nite thermal conductivity in Id lattices |^ and directed 
current in a Id system on nonequilibrium phase tran- 
sition , and on violation of the fluctuation-dissipation 
equality among others. However many theoretical 
works done on nonequilibrium phenomena, to different 
extent, assume certain mathematical forms as the start- 
ing point of their investigations, which, though making 
the problems at hand mathematically attackable, may 
miss generic properties of real physical systems because 
of the limitation of the forms being used. On the other 
hand, with the development of modern computer technol- 
ogy, numerical experiments have become a very powerful 
tool. Here we carry out our investigation by applying 
extensive numerical studies on a simple model system, 
and try to form an understanding based on these results, 
without introducing any ad hoc mathematical functional 
forms to characterize the system. 

We study a two-dimensional system of identical non- 
interacting hard disks confined in a thin pipe (Fig. |^). 
The width of the pipe is set so that two disks cannot 
pass each other. Thus the motion of the disks is two 
dimensional, while at the same time their sequence is 
preserved. The two side walls are periodic and the two 
end walls are energy sources kept at the same tempera- 



ture. For more details, see |9|,[l0|]. We use the simplest 
collision model — after a collision between two disks, the 
normal relative velocity changes sign, and decreases by 
a factor of the restitution coefficient r, with < r < 1. 
In the collision, the other components of the velocities 
are unchanged. Thus the total momentum is conserved 
while a portion of kinetic energy is lost when r < 1 . Sim- 
ulations are for 80 disks and the average spacing between 
two neighboring disks are much larger than their radii. 
Statistical analysis is done for the steady state. Different 
versions of one-dimensional diffusion models are used in 
However, they rely on preset diffusion equations 
as their starting point, while we are attempting to find 
out the physical rules which will lead to those equations. 

Even with this simple model, we nevertheless limit 
ourselves to situations close to equilibrium. Extensive 
works have been done on such systems where linear devi- 
ations from equilibrium are generally observed and stud- 
ied P-^,pd| . In our system, the details of the setting 
are so clearly defined that we may attempt to obtain a 
thorough understanding of the linearity, which then may 
point to the direction for further study on situations far- 
ther away from equilibrium. 

One of the cornerstones of classical kinetics theories, 
Liouville's theorem states that the phase space volume is 
a conserved quantity pd]-p^ . The importance of under- 
lying geometry of the phase space is especially stressed in 
iQ. It is long recognized that energy dissipation makes 
phase space tend to shrink We would like to investi- 
gate the effects of this deviation from Liouville's theorem 
on the dynamics. The set up of the system is to create a 
very small, yet uniform dissipation and look for its first 
order effects. To stay away from disturbing the internal 
dynamics, we put energy sources only at the end walls. 
However, in this way we cannot have energy dissipation 
without energy diffusion. There are two processes going 
on in the system. First, energy dissipation due to inelas- 
tic collisions, which to the first order is uniform. Second, 
energy diffusion, meaning energy conduction from the 
end walls into the system, which at a certain location is 
proportional to its distance to the center — at the center 
there is no net diffusion in either direction due to sym- 
metry. Energy conduction also exists in a similar system 



with elastic collisions between disks but two end walls at 
different temperature, which we will treat later. The co- 
existence of diffusion and dissipation is a general feature 
in dissipative systems. We will demonstrate that they 
can be separated due to different spatial parities they 
exhibit in our system. 
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FIG. 1. A snapshot of the thin pipe system. The periodic 
side walls are indicated by dashed lines. The coordinate sys- 
tem is set up so that the a::-axis is along the pipe. 
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We concentrate on correlations among ai-components, 
i.e. along the pipe, of velocities as an indicator of internal 
orders. For a system at equilibrium, even when there are 
interactions among its constituents, their velocities obey 
Gaussian distribution and thus are uncorrelated. The in- 
fluence of a particle's motion on other particles' velocities 
is an indicator of order, and we quantify this order with 
velocity correlations. Specifically, among other things, 
we look at {ui\vj), meaning the conditional time average 
of the x-component of velocity of the «-th disk when Vj 
is given, where Vj = Uj+i — Uj. Of course, (ui) by itself 
vanishes because in the statistical steady state there is 
no average drift in disks' positions. 

A somewhat more simplified quantity is the collision 



average |10f| , denoted by 



which is the average of Ui 



when the j-th disk is colliding with the {j + l)-st disk. 
For a system in equilibrium, this quantity vanishes when 
i is neither j nor j + 1. However, in our system, it does 
not vanish, but exhibits regular and interesting patterns. 
We can plot ( function of i for several different 

j's — the resulting curves are confusing. But if we plot the 
odd and even parity components of these curves, they are 
easily identified as dissipation and diffusion, respectively. 

The odd parity component, i.e., {ui — U2j+i-i)j is 
for dissipation — they overlap for different value of j's 
(Fig. H), a clear signature for uniform dissipation. Appar- 
ently, dissipation due to collision between two neighbor- 
ing disks has effects on the motions, not just of several 
disks close to them, but of all the disks in the system. 
The even parity component, i.e., {ui + U2j+i-i)j is for 
diffusion — they are proportional to the distance between 
j and the center j — 40 (Fig. H), a signature that they are 
caused by energy conduction which is also proportional 
to this distance. 

For a further test, we run simulations for a system with 
similar set up, but with elastic collisions and the two end 
walls are kept at different temperatures so that no dissi- 
pation, only diffusion occurs. The odd parity component 
of the curves vanishes. And {ui^j)j is independent of 
j, corresponding to uniform diffusion. Also {ui)j can be 
rescaled to overlap with the even parity components from 
the dissipative system (Fig. ^). 
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FIG. 2. The odd parity component of {ui)j. * is for j = 40, 
o for j = 35 and -I- for j — 25. The horizontal axis is for i— j 
rescaled so that the closest boundary to j is shown at 40. Only 
the positive half of i — j is shown. The system has 80 disks, 
r = 0.995 and the two end walls are set to the same tempera- 
ture. The vertical axis is without rescaling, demonstrating a 
uniform correlation structure from uniform dissipation. The 
dashed line is the corresponding rescaled curve for a system 
of 40 disks and r = 0.99. 

Heuristically, we can explain the different spatial pari- 
ties by pointing out that dissipation consists of time irre- 
versible microscopic processes while diffusion consists of 
only time reversible microscopic processes. In the system 
with no dissipation but only constant diffusion, when the 
j-th disk is colliding with the (j + l)-st disk, let us con- 
sider the following two operations: time reverse and a 
spatial parity operation that reverses the direction of x- 
axis and interchanges i-th disk with (2j + 1 — i)-th disk. 
These two operations leave the energy flux unchanged 
and connect two equally probable microscopic states. Be- 
cause Ui in one state is the same as U2j+i-i in the other 
state, we see the correlation curve from diffusion has even 
spatial parity. The odd parity for the correlation curve 
from dissipation is intuitively reasonable, given dissipa- 
tion causes the shrinkage of phase space volume. 

The correlation curves from diffusion and dissipation 
have following properties. They are independent of the 
width of the pipe. The diffusion curves are not sensitive 
as to where the boundaries are — curves for different colli- 
sion number j's can overlap though the distances from j's 
to the boundaries are different. However, they are not for 
local structures. The inset of Fig. || shows the overlap- 
ping between diffusion curves for systems with 80 disks 
(solid line) and 40 disks (dashed line), only after suit- 
able rescaling of i — j, which suggests they are roughly 
functions of (i — j)/N. This is true also for dissipative 
curves, which only overlap with one another after rescal- 
ing of i — j (Fig. g) . These observations suggest that the 
structures from dissipation and diffusion are global. 
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FIG. 3. The even parity component of {ui)j. * is for j = 35, 
o for j = 30 and + for j = 25. The horizontal axis is for i~ j 
without rescaling. Only the positive half of i — j is shown. For 
j = 35, the vertical coordination plotted is the real number 
divided by 5 for the energy flux at j — 35 supplies dissipation 
caused by 5 disks; for j — 30, it is the real number divided 
by 10; for j = 25, it is the real number divided by 15. These 
three curves overlap, showing that the even mode is propor- 
tional to the distance between j and the center. The solid 
line is the corresponding rescaled curve for an elastic system 
with only energy diffusion; its i — j value is rescaled. The 
inset is a similar plot, with rescaled i — j as the horizontal 
axis, also shown as the dashed line is a corresponding curve 
from a elastic system with 40 disks. 

In the rest part of this letter, we describe our detailed 
investigation of correlation curves for diffusion. For this 
purpose, the simulations are for systems with elastic col- 
lisions and different temperatures at two ends. We look 
at the conditional time average of Ui when Vj and Xj are 
given, where Xj is the spacing between j-th and (j-l-l)-st 
disks, and find (Fig. |), 

{u,\Vj,Xj) = C,-j - + f{Xj - {Xj))] , (1) 

where / is a constant. This suggests that the effects of 
Vj on Ui are independent of those from Xj and we can 
concentrate on just the velocities: 

=c,_,(i;f (2) 

It is not difficult to understand this relation for i — j- 
Let us consider {uj -I- Wj+i as a function of Vj. It can- 
not be a constant because = (u'j^i—u'j) < 
while (vj) = 0. It cannot have a linear term propor- 
tional to Vj either because when j-th disk collides with 
(j -|- l)-st disk, Vj changes sign while Uj + is un- 
changed due to momentum conservation. Then the sim- 
plest form {uj + can have is a constant times 
v'j — (v'j), as in (||). Though conforming to quite general 
arguments, the above observation is in contrast to the 



"molecular chaos" assumption of the Boltzmann equa- 
tion, the assumption that makes the equation time irre- 
versible. 




FIG. 4. {ui2\v40,X4o). Though noisy, the surface plot 
shows a linear dependence on X40 and a quadratic dependence 
on 1140. 



It is counterintuitive that does not vanish when 
i is neither j nor j -I- 1. The evidence shown in Fig. ^ 
implies a global structure yet depending on local inho- 
mogeneity. We propose an approach based on geometric 
distortion of the phase space. Let us think the phase 
space as a Riemannian manifold with a metric g for the 
usual {p, q) coordinate system . The statistical prob- 
ability density p is a scalar. But the probability of the 
system in the neighborhood of a point is proportional to 
p^/g■ One approach would be to view the phase space 
as a simple Euclidean manifold with a trivial then p 
has to account for the interesting but puzzling correlation 
structure. Another approach, which we prefer, is to set 
p to a constant, and view the correlation pattern as an 
indication of the distortion of the phase space and thus 
a non-constant g. 

In this approach, because the distortion of the phase 
space, the variables u^'s are no longer statistical indepen- 
dent and the corresponding coordinate system no longer 
orthogonal. Let us assume a new set of independent vari- 
ables u'^s deviates from the original set to the first order: 

u'^ = u, + eF{T), (3) 

where e is proportional to the small inhomogeneity, and 
F denotes the set of Ui's. Because {Ui\vj) = by defi- 
nition, we have {ui\vj) — ~e{F{T)\vj)o. The subscript 
in the last expression indicates that only the zeroth 
order is needed, which can be extracted from numerical 
data. From the somewhat noisy data, we find (H) takes 
the form. 



N-1 N-2 

where C is a constant that keeps (u^) = 0. Let us for the 
moment ignore boundaries and write the coefficients as 
a(i — j) and b(i — j) respectively. b{i—j) can be obtained 
by noticing: 

{ui\vj^i,Vj) - {u^\vj^i,-Vj) = -2b{i - j)vj^iVj. (5) 

Then in turn a{i — j) can be obtained by substituting 
b(i ~ j) into 



-a{i - j) - \a{i - {j - 1)) - ^a{i - (j + 1)) 



(6) 



We use numerical data from (ui|w4o, f4i) to extract values 
of a(i — j) and b(i — j) and plot the result in Fig. |[ 
We can assume the metric under the coordinate system 
diagonal. Then by the relation between two coordi- 
nate systems we can obtain the metric g under the 
original coordinate system of u^'s. Because the metric 
g is connected to the strain tensor we see how the 
correlations indicate the distortion of the phase space. 
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FIG. 5. Coefficients a{i — j) and b{i — j) in (|^. Solid line 
is for a and dashed line for b. The first data point of a is not 
shown, which at about —0.01. The dotted line is a — 2.56 and 
is close to 0. 

Because we study the behavior of the system very close 
to equilibrium, we expect a smooth distribution function 
in the phase space, agreeing with the approach of Tuck- 
erman et al. [|l4| . However, when the inhomogeneity is 
stronger [ p^ . there is likely to be a transition in the dy- 
namics and the distribution function becomes fractal as 
suggested by Hoover, et al. 



In this letter, we present numerical results of velocity 
correlations in a simple model of nonequilibrium system 
and propose to interpret the observation as an indica- 
tion of geometric distortion of the phase space. However, 
most theoretical reasonings in the letter are heuristic ar- 
guments. We did not attempt to formulate a fundamen- 
tal theoretical framework because we feel further data 
collecting, and more importantly, different, more general 
setting of the system should be tested before a better un- 
derstanding can be reached. A simple extension would be 
the testing of a similar system, but with two dimensional 
close-packed disks. 
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